Dissipative soliton excitability induced by spatial inhomogeneities and drift 
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We show that excitabiUty is generic in systems displaying dissipative sohtons when spatial inho- 
mogeneities and drift are present. Thus, dissipative solitons in systems which do not have oscillatory 
states, such as the prototypical Swift-Hohenberg equation, display oscillations and Type I and II 
excitability when adding inhomogeneities and drift to the system. This rich dynamical behavior 
arises from the interplay between the pinning to the inhomogeneity and the pulling of the drift. 
The scenario presented here provides a general theoretical understanding of oscillatory regimes of 
dissipative solitons reported in semiconductor microresonators. Our results open also the possibility 
to observe this phenomenon in a wide variety of physical systems. 

PACS numbers: 05.45.Yv, 42.65.Sf, 89.75.Fb, 02.30.Jr 



Spatially extended ^sterns display a large variety of 
emergent behaviors including coherent structures. 
Particularly interesting is the case of dissipative solitons 
(DS) exponentially localized structures in dissipative 
systems driven out of equilibrium, as they can behave like 
discrete objects in continuous systems. DS can display a 
variety of dynamical regimes such as periodic oscillations 
, chaos [1, 0] J or excitability ^ . DS emerge from a 
balance between nonlinearity and spatial coupling, and 
driving and dissipation. They are unique once the sys- 
tem parameters are fixed, and they are different from 
the well known conservative solitons that appear as one- 
parameter families. 

In optical cavities, DS (also known as cavity solitons) 
have been proposed as bits for all-optical memories 0- 
[isj . due to their spatial localization and bistable coexis- 
tence with the fundamental solution. In Ref. [81 it was 
reported that DS may exhibit excitable behavior. A sys- 
tem is said to be excitable if perturbations below a cer- 
tain threshold decay exponentially, while perturbations 
above this threshold induce a large response before go- 
ing back to the resting state. Excitability is found for 
parameters close to those where a limit cycle disappears 
[14| . Excitability mediated by DS is different from the 
well known dynamics of excitable media, whose behavior 
stems from the (local) excitability present in the system 
without spatial degrees of freedom. Excitability of DS is 
an emergent behavior, arising through the spatial inter- 
action and not present locally. Moreover, the interaction 
between different excitable DS can be used to build all- 
optical logical gates [15j. 

In real systems, however, typically solitons are static, 
so oscillatory or excitable DS are far from being generic. 
In this work we present a mechanism that generically 
induces dynamical regimes, such as oscillations and ex- 
citable behavior, in which the structure of the DS is pre- 
served. The mechanism relies on the interplay between 



spatial inhomogeneities and drift, and therefore can be 
implemented under very general conditions. Inhomo- 
geneities, or defects, are unavoidable in any experimen- 
tal setup, and drift is also often present in many optical, 
fluid and chemical systems, due to misalignments of the 
mirrors 16|, Il7| , nonlinear crystal birefringence 18| , or 
parameter gradients [31 ) in the first case, and due to the 
flow of a fluid in the others 2^, 21|. Roughly speaking, 
the presence of inhomogeneities and drift introduce two 
competing effects. On the one hand, an inhomogeneity 
pins a DS at a fixed position and, on the other, the drift 
tries to pull it out. If the drift overcomes the pinning, 
DS solitons are released from the inhomogeneity 3 22 1. 
Thus, depending on the height of the spatial inhomogene- 
ity and the strength of the drift, we observe three main 
dynamical regimes: i) stationary (pinned) solutions, ii) 
oscillatory regimes, where DS source continuously from 
the inhomogeneity, and iii) excitability, where a pertur- 
bation may trigger a single DS that is driven away from 
the inhomogeneity. This scenario does not depend on 
the details of the system nor on the spatial dimension- 
ality. Besides, it provides a solid theoretical framework 
to explain the dynamics of DS in systems with drift, as 
observed for instance in semiconductor microresonators 
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To analyze this scenario we start with the prototypi- 
cal Swift-Hohenberg equation (SHE), a generic amplitude 
equation describing pattern formation in a large variety 
of systems The SHE is variational, thus it can 

not have oscillatory regimes. Adding drift and a spatial 
inhomogeneity we have: 
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Here u{x,t) is a real field, a and g real parameters, and 
c is the velocit y o f the drift. The SHE displays DS for 
a > V27/38g QUil, 



so we take a = 1.2 and 17=1. 
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FIG. 1. (Color online) Bifurcation diagram showing the max- 
imum of the steady state, Umax, as a function of h for c = 0. 
In black solutions pinned at their center and in gray (blue) 
DS pinned at the first tail minima. Solid (dashed) lines indi- 
cate stable (unstable) solutions. The insets show the profiles 
of the solution and the spatial inhomogeneity. 



We fix fcp — 0.5. We use a supergaussian gain profile 
r{x) = ro — 1 -l-exp{ — [{x ~ xq) /e]^^} to model a finite sys- 
tem of width 2e so that DS disappear at the boundaries. 
We fix ro = —0.2, roughly in the middle of the subcrit- 
ical region where DS exist, and, except where otherwise 
indicated, e = 94.0842. Small changes of the parameter 
values do not substantially modify the results. A single 
spatial inhomogeneity located at the center is introduced 
adding a Gaussian profile b{x) = /iexp{— [{x — X())/(j]'^} 
with height h and half- width a [2^ . This is motivated by 
the addressing beams used in optical systems to control 
DS [13, E^l- We take a = 2.045, roughly half the width 
of a DS. Similar results are found for other values of a 
provided the width is small enough to avoid trapping two 
DS. We have also checked that similar results are found 
for a defect in the gain profile r. We take h and c as 
control parameters [28]. 

For c = we observe three different steady state solu- 
tions which are the main attractors of the dynamics (Fig. 
[T]): i) the fundamental solution, a low bump correspond- 
ing to the deformation of the homogeneous solution, ii) 
a high amplitude DS pinned at its center, and iii) a DS 
pinned at the first oscillation of its tail. Decreasing h, at 
/i = 0, a transcritical bifurcation occurs in which branch 
ii) becomes unstable while iii) is stabilized. Physically, at 
h = the defect goes from being a bump to a hole. DS 
tend to sit at the inhomogeneity maximum, thus a DS 
centered at the hole becomes unstable and shifts its po- 
sition until the hole coincides with the first minimum of 
its tail [2^ . At /i = there is also a crossing of unstable 
middle-branch DS. In what follows we focus on h > 0. 

This scenario changes when drift is introduced, c ^ 0: 
parity symmetry, x —x, is broken making the trans- 
critical bifurcations imperfect. This implies a rearrange- 
ment of branches, that for low values of c leads to the 
snake- like branch shown in Fig. ^jji.). Saddle- node bifur- 
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FIG. 2. The same as in Fig. [T] for a) c = 0.05, b) c = 0.12 
and c) c = 0.4 (see main text). In c) the crosses indicate 
the maximum and minimum value of the oscillatory DS at a 
given spatial location. The grey dotted line is a sketch of the 
unstable cycle for illustrative purposes. 
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FIG. 3. a) (Color online) Two-parameter (c vs K) phase dia- 
gram of the system (see text). Panels b) and c) show a zoom 
close to the cusps and Takens-Bogdanov points. 



cations SNi and SN3 where already present for c = 0. A 
Saddle-Node on the Invariant Circle (SNIC) reconnects 
branches ii) and iii), while the saddle-node SN2 arises 
from the middle branches transcritical bifurcation. In- 
creasing c the branch stretches (cf. Fig. [2j)) and SNi 
coalesces with SN2 at the cusp bifurcation Ci as shown 
in Fig. [3] which displays the general scenario in parameter 
space. Close to Ci there is a Takens-Bogdanov point TBi 
that unfolds a Hopf bifurcation line H+ [30,] (Fig. [SJj) . As 
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FIG. 5. Excitable excursion of the fundamental solution in 
region I, close to FC (c = 0.6, h = 0.085), triggered varying 
h an amount Ah = 0.0325 during At = 20 time units. The 
dot-dashed line shows the gain profile r{x). 



FIG. 4. a) (Color online) Snapshot of a train of solitons in 
region II of Fig. [3]for c = 0.11 and h = 0.076. The solid grey 
(red) and the dot-dashed lines show b{x) and r{x) respectively. 
Here e = 163.625. b) Time evolution of the field u at xo- c) 
Field profile at the times shown with symbols in b) showing 
how a new DS develop and detach. 
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c keeps increasing, the SNIC line encounters a Saddle- 
Node Separatrix Loop (SNSL) codimension 2 bifurcation 
from which a saddle-node SN4 and a saddle-loop SL bi- 
furcation lines unfold (Fig. [3]:). SN4 soon coalesces with 
SN3 at the cusp C2 while SL ends at a Takens-Bogdanov 
point, TB2, which also unfolds a Hopf line H~ tangent to 
the SL line. Finally for larger values of c there is a single 
monotonous branch of steady state solutions (Fig. [2t). 
The Hopf line H"*" is subcritical and is accompanied by 
a fold of cycles (FC) from which a stable cycle emerges. 
For h large enough the defect is above the threshold to 
switch on a DS. Thus this limit cycle corresponds to the 
periodic creation of DS at the inhomogeneity that are 
then drifted away (Fig. |4};), generating a train of soli- 
tons that, in our case, disappear at the boundary of the 
system (Fig. 0^) . Such regime was observed experimen- 
tally in semiconductor microresonators , where it was 
called a "soliton tap" . This cycle is stable all the way to 
the supercritical Hopf H~ , where h becomes large enough 
to prevent the DS being advected by the drift, and the 
cycle ends. 

In parameter space (Fig. [3]) the fundamental solution 
is stable in regions I and V, pinned DS are stable in HI 
and IV and stable trains of DS exist in II and V. The sce- 
nario shown in Figs.Et" andO:, much richer than in Q, 
is characteristic of systems displaying relaxation oscilla- 
tions described, for instance, by the Van der Pol equation 
(see Fig. 2.1.2 in 3l|). These systems are characterized 
by two different time scales, a slow and a fast one. Here 
the fast time scale corresponds to the switching time of a 
DS on the inhomogeneity, while the slow one is the time 
the drift takes to detach a DS once it is formed [22]. The 
switch-on time is basically independent on the drift, while 
the detach time has a strong dependence on it. These two 
time scales can be clearly identified in a time trace of the 
field at Xo (Fig. Hb) . Excitability can be found in regions 



FIG. 6. Oscillation period T in region II as a function of h 
close to a) the FC (c = 0.6) and b) the SNIC (c = 0.12). The 
inset in b) shows the scaling of the divergence of the period. 



I, HI and IV when applying a perturbation for a short 
time that changes either c or h. If the system is brought 
monetarily into the oscillatory region II then an excitable 
excursion consisting on the exploration of a loop of the 
cycle before returning to the original steady state is trig- 
gered. An alternative way to trigger an excitable excur- 
sion is perturbing the state of the system rather than a 
parameter 14]. Transient parameter changes are usually 
easier to implement experimentally, thus in what follows 
we consider this kind of perturbations. 

In region I a superthreshold perturbation of the fun- 
damental solution grows to generate a DS, that is then 
advected away, setting the defect back to the resting state 
(Fig. [5]). As the transition from the stationary state to 
the oscillatory one is mediated by a subcritical Hopf bi- 
furcation, the excitability is of Type II ^14] • A clear sig- 
nature of this type of excitability is that the period of 
the oscillations remains practically constant as one ap- 
proaches the threshold from the oscillatory side, namely 
as one approaches FC coming from region II, as shown in 
Fig. ISK- As typical for Type II excitability the threshold 
is rather a quasi-threshold. If one applies a perturbation 
that crosses FC but not H+, the system will be stacked 
in region V (bistable) . The closer the parameters are to 
FC, and the closer is FC to H+, the smaller is the per- 
turbation required to trigger an excitable excursion [32] . 

Different dynamics is found in region HI, where h is 
large enough that, beyond switching on a DS if there is 
none, it also pins it once it is formed. In this case DS un- 
dergo a supercritical Hopf bifurcation when crossing H^ . 
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FIG. 7. The same as in FigObut close to H , in region III. 
Here c = 0.4, h = 0.18, and Ah = -0.045. 



Close to threshold, DS exhibit small oscillations, but de- 
creasing h or increasing c just a little further DS start to 
source from the inhomogeneity, forming a sort of canard 
in phase space [s^ , leading to the train of DS found in 
region II. This is again a mechanism leading to Type II 
excitability as shown in Fig. [T] At difference with the 
previous case, here the resting state is a high amplitude 
DS (see first panel), and the excitable excursion consists 
of the DS leaving the inhomogeneity and a new one being 
formed. 

There is a third mechanism, associated to the SNIC 
line separating regions II and IV. A detailed analysis of 
DS excitability mediated by a SNIC can be found in Ref. 
[26| . For parameters in region IV a supra-treshold per- 
turbation that brings the system into region II triggers 
the unpinning of a DS leading to an excitable excursion. 
The initial and the final state is a pinned DS and the ob- 
served behavior is very similar to the one shown in Fig.m 
In this case excitability is of Type I, whose signature is a 
divergence of the period. As shown in Fig. IBb) the period 
diverges with a power law of exponent —1/2 approaching 
the SNIC line from region II. A divergence in the period 
has been reported in semiconductor microresonators .221] , 
induced also through spatial inhomogeneities and drift. 
Additionally, the SL line in Fig.|3t) leads to another Type 
I excitable region (discussed at length in Q), but at dif- 
ference with the SNIC bifurcation line, that extends for 
a wide parameter range, the SL occurs in a very narrow 
region (Fig. [3t), and therefore it may be more difiicult 
to find experimentally. 

In summary, we have shown that the competition be- 
tween the pinning of DS to a spatial inhomogeneity and 
the pulling generated by a drift leads to a complex be- 
havior leading to oscillations and to Type I and II ex- 
citability through several mechanisms. To clearly illus- 
trate this we make use of the SHE, for which there is a 
Lyapunov potential thus the system always evolves 
to steady states that minimize the potential. When drift 
and boundary conditions are added, the SHE no longer 
has a Lyapunov potential 
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has been recently used to describe DS observed experi- 
mentally [Hill. Our results do not depend on the mi- 
croscopic details of the system, or the number of spatial 
dimensions, or the nature of the spatial inhomogeneity, 
but on general emergent properties of DS, providing a 
theoretical framework to explain the dynamics of DS in 
presence of spatial inhomogeneities and drift in a very 
broad class of extended systems. The theory presented 
here gives an explanation to experimental observations 
[i^ l and opens up the possibility to observe this phe- 
nomenon in a variety of other experimental setups. 
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behavior can arise. In order to confirm the generality of 
our results, we have checked that the scenario is qualita- 
tively the same in a completely different model displaying 
DS [351] , namely the Lugiato-Lefever equation i36] , which 
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